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A simple description of superfluid hydrodynamics in the inner crust of a neutron star is 
given. Particular attention is paid to the effect of the lattice of nuclei on the properties of 
the superfluid neutrons, and the effects of entrainment, the fact that some fraction of the 
neutrons are locked to the motion of the protons in nuclei. 

§1. Introduction 

The outer layers of a neutron star at densities below nuclear saturation density 
make up only a small part of the star, but they are important because processes 
that occur there are more accessible to observation than are those in the inner parts 
of stars. Among phenomena for which the physics of the outer layers of the star is 
crucial one may mention the thermal evolution of neutron stars. This includes the 
recent suggestion that heat transport by superfluid neutrons in the crust could be a 
significant contributor to thermal transport!- 1 -' Another example of the importance 
of the crust is the proposal that quasiperiodic oscillations seen in x-ray afterglows in 
giant flares on neutron stars may be attributed to excitation of torsional modes of 
oscillation of the crustP~® j n both these contexts, modes of oscillation in the crust 
play a key role, and the purpose of this article is to present a simple approach to 
calculating their properties. 

The focus of this article is the inner crust of neutron stars, where superfluid neu- 
trons coexist with a lattice of ions. There is an extensive literature on the physics of 
this region (e.g!^ 1 ). However it has been only recently pointed out that an important 
role is played by "entrainment" of neutrons by the protons in nucleiP 1 Such mutual 
entrainment effects have been studied for a long time in the liquid core of neutron 
stars (see e.g. 71 ). In this article, we begin by describing how to formulate the 
theory of superfluidity when there is a periodic lattice of ions. We then present the 
hydrodynamic equations and show that, as a consequence of Galilean invariance, an 
appropriately defined superfluid density characterizes both the effective mass of nu- 
clei and the strength of a vector coupling between lattice vibrations and sound waves 
in the neutron superfluid. We derive the dispersion relation for the long-wavelength 
collective modes, which are generally a superposition of lattice vibrations and sound 
waves in the neutron superfluid. Finally, we comment on the physical meaning of 
the neutron superfluid density. 
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§2. Neutron superfluidity in the presence of a lattice of nuclei 

In a Fermi superfluid, the pairing amplitude at position r is given by (ip^(r)ip^(r)) , 
where Vv( r ) is the annihilation field operator for a neutron with spin a at position 
r, and (. . .) denotes a quantum-mechanical or, at nonzero temperature, a thermal 
average. For a superfluid at rest, the phase 2(p(r) of the pairing amplitude is in- 
dependent of position, while for a superfluid in motion it depends on space. In a 
translationally invariant system with paired neutrons, such as neutron matter in the 
outer core of a neutron star, the phase of the pairing amplitude varies smoothly in 
space, while in the presence of a crystal lattice, the phase has, in addition, oscillatory 
contributions with a period equal to the lattice spacing. The question is know what 
quantity should one use to describe motions of the superfluid on length scales large 
compared with the lattice spacing. 

Just as for single-electron wave functions in a periodic lattice, one may char- 
acterize states of the neutron superfluid in the crust by a wave vector k such that 
the change of the phase y(r) when the spatial coordinate is displaced by a lattice 
vector R is given by k • R. Alternatively, one may work in terms of a coarse-grained 
average phase, 0(r), which is the average of <p over a region in the vicinity of r with 
dimensions large compared with the lattice spacing but small compared with other 
length scales in the problem, in which case k = V</>. Similar physics occurs for atomic 
Bose-Einstein condensates in artificial crystal lattices created by standing-wave laser 
beams, so-called optical lattices.® 

§3. Hydrodynamic equations 

We now derive equations of motion for the superfluid and the crystal lattice 
in situations when the scale of spatial variations is large compared with the lattice 
spacing. In addition, we shall require that the scale of spatial variations is large 
compared with the size of the "neutron pair wave function", which is given by the 
superfluid coherence length, £ ~ hvFn/A n , where vp n is the neutron Fermi velocity 
and A n the neutron superfluid energy gap. We shall treat the nucleons as nonrela- 
tivistic, which is a good first approximation in neutron-star crusts. For simplicity, 
we shall neglect the effects of the gravitational potential. Since this varies relatively 
little over the crust, it is a reasonable approximation to treat space as flat. 

Hydrodynamic equations are an expression of conservation laws for momentum 
and the numbers of neutrons and protons. Our treatment follows closely the stan- 
dard discussion of the two-fluid model for liquid helium-4P For simplicity, we shall 
assume here that the temperature is zero and we shall neglect the external potential 
and inhomogeneity of the equilibrium density distribution. We shall consider distur- 
bances whose length scale is large compared with the screening length for electrons, 
in which case matter is electrically neutral locally, and the number density of elec- 
trons is equal to that of the protons. Two of the variables with which we choose 
to work are two densities - the total neutron density n n and the proton density 
n p . This choice is particularly convenient, since these are the quantities that enter 
the conservation laws for particle number. Alternatively, one could have chosen use 
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appropriately defined densities of neutrons associated with nuclei and with the "su- 
perfluid between nuclei" P 1 This does not alter the final results, but the formalism is 
particularly straightforward when n n and n p are used. The two other variables we 
employ describe the motion of the constituents. We shall assume that the ions are 
sufficiently far apart that hopping of protons from one nucleus to another is absent, 
so we may describe the configuration of the ions in terms of the displacement £(r) of 
an ion from its equilibrium position. Thus the velocity of the ions is given by v p = 
For the superfluid, we characterize its motion by the wave vector k, or equivalently 
the "neutron superfluid velocity" , which we define as^ 1 

Kk 

v n = — , 3-1 

m 

where m is the nucleon mass (we neglect the difference between the neutron and pro- 
ton masses). The choice of the mass used in this expression is a matter of convention, 
since the important quantity is really the wave vector. The first two conservation 
laws we consider are those for conservation of the number of neutrons and the num- 
ber of protons. In the case of protons, the continuity equation, when linearized, 
is 

+ n p V • v p = 0. (3-2) 

Throughout, we shall consider small departures of the system from equilibrium, 
and therefore we shall linearize the equations. The full non-linear equations can 
be found in Refs. EJ™ We shall assume that the atomic number Z of nuclei is a 
constant in the region under consideration, and therefore the proton density is given 
by n p = Znj, where nj is the density of nuclei. To write a conservation law for the 
neutron number, we need an expression for the current density of neutrons. Consider 
first the case where the ions are at rest, and we denote quantities in this frame by a 
prime. Phenomenologically, for small superfluid velocities one may write the current 
density of neutrons as 

j'n = < v n, (3-3) 

where n s n may be regarded as the density of superfluid neutrons. Now let us consider 
the current in a frame moving with velocity — v p . By Galilean invariance, the current 
in the new frame is given by 

in = n s n V n + n n v p = n s n v n + n b n v p . (3-4) 

Here v„ = v' n + v p is the neutron superfluid velocity in the new frame. This result 
follows from the fact that under the Galilean transformation, every single nucleon 
state acquires an extra phase factor exp(imv p -ri/h), r-j being the nucleon coordinate, 
and therefore the phase (/?(r) is increased by mv p ■ r/h . The quantity n n = n n — n s n 
may be regarded as the density of neutrons bound to nuclei. The continuity equation 
for neutrons is therefore 
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The next equation is that for conservation of momentum, which in the absence of 
dissipation reads 

| = f, (3-6) 

where / is the force per unit volume. The momentum density may be calculated in 
a similar way to that in which we calculated the current above. If effects of relativity 
on the motion of the neutrons are neglected, the momentum density is given by the 
current density, times the nucleon rest mass, and thus, in the frame in which the 
protons are at rest, 

9n = m < v n- (3'7) 

By making a Galilean transformation to a frame moving with velocity — v p , one 
finds that the momentum density is changed by an amount pv p , and therefore the 
momentum density is given by 

9 = /0v p + mn s n v' n = pfv p + mn s n v n . (3-8) 

Here p is the total mass density, the energy density divided by c 2 , and p^ = p — mn s n . 
If the electron contribution to p is neglected, p eS = m(n p + n^), and the momentum 
density is the nucleon current density times m: 

g = mn s n v n + m(n p + n b n )w p . (3-9) 

The condition for momentum conservation, (j3-6|) . is therefore 

m<-^ + m(n p + n b n )-^ = f. (3-10) 

When the rigidity of the crystal lattice is neglected, the system may be treated as 
coupled fluids and the force is given by 

f = -Vp, (3-11) 

where p is the pressure. At zero temperature, changes in pressure are related to 
changes in the neutron and proton chemical potentials by the Gibbs-Duhem relation 

dp = n n dp n + n p dp p , (3-12) 

and the condition for momentum conservation becomes 

mn n-gf + m ( n P + n n">~gf + n n^ + n p V (J, p = 0. (3-13) 

When the protons form a rigid lattice of nuclei, it is necessary to include the 
effects of elasticity. For an isotropic solid, the elastic contribution to the energy is 
given to second order in the strain tensor 

1 / dui dun 



- - 



2 \dxj dxi 



+ id- (3-14) 
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E el3Stic = -Kui + S[ Uij -^-^\ . (3-15) 



where K = d 2 E/du 2 i \ rin and S is the shear modulus. In a solid with lower symmetry, 
the expression for the elastic energy has more terms. For crustal material with 
roughly spherical nuclei, the crystal symmetry is expected to be body-centered cubic. 
However, it seems unlikely that the crust behaves as a perfect crystal, and one would 
expect it to be polycrystalline. On length scales large compared with the size of a 
domain, the solid will behave as an isotropic solid. For phases with liquid-crystal-like 
phases, the elastic properties are even more complicated, and the effects of imperfect 
ordering of nuclei has yet to be investigated.^ 1 

In a neutron star crust, there is also a contribution due to coupling of the 
neutrons to the strain field, and the contribution linear in the strain and linear in 
the changes in the neutron density has the form 

^lattice = —L5n n ua, (3- 16) 

which is the only term consistent with rotational invariance. The trace uu of the 
strain tensor is related to changes in the proton density by the equation 

u n = - 6 -^. (3-17) 

n p 

Thus one can see that 

L _ d 2 E OfXn ( - 348 x 

dn n dua p dn p 

Similarly, one can show that K = n 2 dfi p /dn p . For an elastic solid, the derivatives 
with respect to n p are to be evaluated with the shear strain held fixed. For a 
solid consisting of a single component, the quantity K would be the bulk modulus, 
whereas for the system under consideration here, it is the contribution to the bulk 
modulus when the density of neutrons is held fixed. The total bulk modulus is 
given by K = —Vdp/dV, where V is the volume of the system and the derivative 
is to be evaluated keeping the numbers of neutrons and protons fixed, and therefore 
K = K + 2n n L + r? n d[i n jdn n . Quite generally, the quantity d/ii/drij represents 
an effective interaction between densities of components i and j. Thus the effective 
proton-proton interaction is K /n 2 and the neutron-proton one L/n p . 

With the effects of elasticity included, the force per unit volume is given by 

fi = -n n ViLi n + KViUjj - LVn n + 2SVjUij, (3-19) 

and the condition for momentum conservation becomes 

mn s n —^- + m(n p + n*)-^ + n„Vj/i n - KViUjj + LVjn n - 2SVjUij = 0. (3-20) 

A fourth condition is needed to close the system of equations. This is the Joseph- 
son relation, which gives the time evolution of the phase of the superfluid neutrons. 
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In a uniform system, the rate of change of the phase is given by 

^ - (3-211 

dt ~ n ' [6 ' 

and provided spatial variations are sufficiently slow, this condition will also apply 
locally. The spatial gradient of this relation gives an equation for the evolution of 
the superfluid velocity which, because of Eq. (|3-ip . is 



<9v 1 

+ -Vfi n = 0. (3-22) 
at m 

If one subtracts n s n m times this equation from Eqs. (|3-13p and (|3-20p one finds 

m(n p + n b n )-^r- + ™«V/i n + n p V^ p = 0, (3-23) 

and 

m(n p + n b n )-^- + n h n V^ n + LV % n n - KV lUjj - 2SV jUij = 0. (3-24) 

Equations (|3-23p and (|3-24p show that nuclei behave as though they had a mass 
density m((n p + n b n ))/ni = Zm(l + n h n /n p ). The quantity n b n thus represents the 
number density of neutrons that should be regarded as being bound to protons, so 
far as the inertia of the nuclei is concerned. 

§4. Normal modes 

As an illustration we consider normal modes of oscillation of the system, which 
are obtained by solving the linearized hydrodynamic equations, treating velocities 
and deviations of the densities from their equilibrium values as small quantities. We 
expand the deviations of the neutron and proton chemical potential by writing for 
the fluid case 

<9/% d\i n 
- — dn n + -— 

on„ on r 



d^n = -^r-dn n + ~^-dn v (4-1) 



and the corresponding equation obtained by interchanging the neutron and proton 
labels. For the case of a solid, the neutron chemical potential is given by jx n = 
dE/dn n \ Uij , and therefore 

d^n = Tr^drin - Lduu. (4-2) 
on n 

The normal modes have the form of plane waves that vary in space and time as 
e i(q.r-u;t)^ w h ere q j s tli e wave vector and cu the frequency. In an isotropic medium, 
the normal modes may be separated into transverse and longitudinal ones. At long 
wavelengths, the normal modes all have a sound-like dispersion relation, with uj = vq, 
v being the mode velocity. The transverse modes are the simpler, since the neutron 
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superfluid does not respond to long- wavelength transverse disturbances, and the 
dispersion relation is given by 

(4-3) 



(n p + n h n )m 



This differs from the usual relation for a one-component solid, v\ = S/p by the 
replacement of the total mass density, p, by (n p + n b n )m, the mass density of protons 
plus neutrons that move with the protons or, equivalently the total mass density less 
the density of superfluid neutrons, p — mn s n . 

The dispersion relation for longitudinal waves has the form 

(mvf - n s n E nn ) (mvf(n p + n b n ) - (K + 45/3)) 

-L 2 < - 2mvfn b n L - mv 2 {n b n ) 2 E nn = 0. (4-4) 

characteristic of two coupled modest One mode is a sound wave in the neutron 
superfluid, which has a velocity (n s n E nn /m) 1 / 2 and the second is a longitudinal lattice 
sound mode with velocity {{K +4S/3)/m{n p +n^ l ) 1 ' 2 . The coupling between the two 
modes has three contributions, the 1? term, which takes into account interactions 
between densities, the term proportional to {n b n ) 2 , which represents a current-current 
interaction, and a mixed term. In the crust of a neutron star at a density just above 
that for neutron drip, the velocity of the sound mode in the neutron superfluid is 
much less than that of the lattice mode, while closer to the inner boundary of the 
crust, the two velocity are similar and there is significant hybridization. Detailed 
calculations will be reported elsewhere. 



§5. The neutron superfluid density 



We may write Eq. (|3-4p for v p = in the form j n = (n^/m)frk. The quantity 
Xjj = n s n /m is the static, long- wavelength current-current correlation function for 
the neutrons, and it thus describes the change in the current density of neutrons 
when the crystal momentum Kk of the superfluid is changed, the ions being station- 
ary. Alternatively, one may write g n = n s n (Kk). The proton mass that enters the 
formalism in §3 is the bare mass, m, since for non-relativistic particles, the current 
density and the momentum density are related by the equation g n = mj n . One may 
write Xjj = n n ff / m n ff > where n^ ff is an "effective density of superfluid neutrons" and 
a "neutron effective mass". If one takes as n° ff the density of neutrons with 
single-particle energies greater than the mean field interaction energy outside nuclei, 
one is led to the conclusion that rrf^ can be as large as ~ 15mP^ However, in the 
hydrodynamic equations, n^ ff and m^ ff will occur in the combination n^ ff /m^ ff , so 
it is irrelevant how one distributes changes in \jj between a change in the effective 
mass and a change in the effective density of superfluid neutrons. The quantity n s n 
encapsulates all information about the superfluidity of neutrons that is needed for 
a hydrodynamic description: it includes effects due to the fact that many bands are 
filled, and that the effective mass of neutrons in conduction bands may be different 
from the bare massP^'^S 1 
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§6. Concluding remarks 



In this paper we have indicated how to describe superfluidity of neutrons in the 
periodic lattice of ions in the inner crust of neutron stars. We have also derived 
hydrodynamic equations to describe collective oscillations. The key ingredients are 
the effective interactions dfj,i/drij,(i,j = n,p) between neutrons and protons, the 
shear elastic modulus S, and the "neutron superfluid density" n s n . The latter quan- 
tity describes how the neutron current density responds to a vector potential, and 
it is not simply related to the neutron density between nuclei. The related quantity 
n\ = n n — n s n determines both the effective mass of nuclei and the vector coupling 
between modes of the superfluid neutrons and lattice phonons. The fact that the 
same quantity occurs in these two contexts is a consequence of Galilean invariance 
of the combined system consisting of the lattice and neutrons. Applications of the 
formalism above will be reported elsewhere. 
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